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This final report is divided into 
three separate parts. The first part 
is a summary of work accomplished ■ : 
during the performance of the con- 
tract. The second part is a report 
on the " Feasibility Study of Laser 
Operations in a Vacuum ", The final 
part is a review of some work in 
progress at PMS relating to the 
"Theoretical Response of Active Scat 
tering Particle Spectrometers", The 
three parts are arranged such that 
they-- -can- -be •'removed - out of context 
and stand alone. 
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PART I 

Summary of Work Accomplished 

The work under this contract was initiated during June, 1972 with 
early work concentrating on the study of laser operation in a vacuum. 

Part II details this study which resulted in a laser design which allowed 
the laser plasma tube and electrodes to be sealed off from the vacuum. 

Thus the whole question of potential electrical and thermal problems 
relating to laser -operation in a vacuum was precluded. The majority of 
the work conducted under this contract was completed on or ahead of schedule 
with delivery a month ahead of schedule. 

The observation of the ASPS in operation brought with it a couple 
of interesting -findings. It *was observed Ihat for small particles the 
scattering signal pulses had some high frequency ripple content. The 
ripple is believed to arise from two sources. One is the Doppler frequency 
created by interference of the forward and backward scattered light which 
can occur when the particle has an axial velocity component. For an 
axial velocity component of one centimeter per second the Doppled velocity 
would be in the tens of kHz range. Frequencies were generally observed in 
the range of 10-100 kHz. However, on occasion frequencies were also 
observed in this range when particles were essentially stationary. This 
may be due to mode switching but we are not at all clear at this time how 
this should arise. Mode switching could not be observed in. the output beam 
intensity distribution nor did the reference beam have within an order 
of magnitude the percentage of modulation that the scattering pulses possessed. 
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The scattering pulses often appeared to have 20% modulation. The fact that 
this ripple frequency is considerably greater than the frequency content of 
the fundamental pulse frequency would allow for it to be separated by filter- 
ing. A correlation study could be performed to determine if the frequency 
or its percentage of modulation were in any way a function of particle size. 
The Doppler effect can be estimated and controlled by varying particle 
trajectories. We believe these observations demand closer scrutiny and study. 
PMS will perform such studies in the future on internal funds if outside 
support cannot be found to further this work. 

The calibration of the two ASPS size ranges involved both theoreti- 
cal and emperical studies. The two system' s theoretical calibration curves 
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of scattering signal versus size closely match an r function. The only 

discrepancy being in the range of 0.1 to 0.3 microns where the curve is 

valid only for high refractive indices. This,.is _no,t .Relieved . problematic 

in that low refractive index materials would be rare in a hard vacuum 

since these are the more volatile species. Thus, no attempt was made to fit 

2 

the pulse height detector response to other than an r function. This also 
gives interchangeability between the pulse height detector circuit boards 
which is the fastest procedure when trouble shooting. 
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PART II 

Feasibility Study of Laser Operations in a Vacuum 
I. INTRODUCTION 

Prior to the design of the ASPS, a study was necessary to determine 
the operational characteristics of open cavity gas lasers under vacuum con- 
ditions. Three questions immediately arise when considering gas laser 
operation under vacuum conditions. First, the laser dissipates a consider- 
able amount of heat which must somehow be dissipated. Second, a gas laser 
is an active interferometer which in the design of such for particle size 
measurements has an extended open cavity. The effect of the change in 
medium — from one atmosphere air pressure to a hard vacuum must be considered 
as having possible effects on laser cavity stability. ..Thirdly, the high 
potentials necessary to start and sustain laser action require protection 
to avoid electrical discharges and corona. Other questions could also be 
considered, such as particle trajectories at velocities under vacuum con- 
ditions. However, these are obviously factors which could not be resolved 
nor influence laser operating characteristics. 

What follows are three separate sections which largely deal with 
theoretical considerations connected with the thermal, interferometric and 
electrical problems, or otherwise influences of laser operations under 
vacuum conditions. The fourth section extracts basic design criteria from 
the previous three sections and summarizes the thoughts behind PMS's 
decision to seal off the laser plasma envelope and its electrodes from the 


vacuum. 
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II. THE HEAT DISSIPATED BY LASER PLASMA TUBES 

Laser plasma tubes of the type used in ASPS systems generally use 
a gas fill mixture which generates about 150 volts/in drop along the discharge 
path length. With an 8" long plasma tube, a voltage drop of 1200 volts would 
be anticipated. A current of Uma would thus generate a power dissipation of 
about 5 watts. 

The Ballast resister is normally placed at the anode of lasers and 

2 

its power dissipation would also cantribute to the total (approximately I R= 

2 

(.004 x 125K = 2 watts). However, we have operated ASPS lasers with 
ballast separations of several feet without difficulty and always contained 
such ballast within the power supply. 

The 5 watts may be radiated in whole, or in part, or conducted away 
by a 'mounting structure. The two extremes to consider for the radiative pro- 
blem were for the laser envelope to be looking at infinity with an effective 
radiative equalibrium temperature of 0° Kelvin or infinity with an effective 
radiative equalibrium temperature of 300° K. The computed radiative equil- 
ibrium temperatures of the laser plasma tube for these two extremes are 376°K 
and 414°K respectively-*-. From this it is obvious that the laser would require 
some conduction to keep the temperature below the softening point of the 
epoxy used to seal the mirrors and Brewster's windows. PMS found that an 
aluminum mount having contact with 60% of the plasma envelope and a conducting 
cross section of 3 cm^ would be required to keep the envelope within 10°C of 
ambient (metal to which tube mount connects is here considered to be at 


L An emissivity of 0.96 was used for glass envelope of 42 cm surface 


area. 
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ambient temperature. It is obvious from the above analysis that the thermal 
problem created by laser operation in a vacuum is an important consideration 


for ASPS system design. 
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III. THE EFFECT OF ATMOSPHERIC PRESSURE ON OPEN CAVITY LASER OPERATION. 

In the analyses of the effects of pressure variations one must con- 
sider two different cavity geometries. The plasma tube best suited for 
ASPS measurements is a hybrid tube with one sealed mirror. The mirror 
configuration is hemispherical (one plane mirror and one spherical mirror) 
which gives a conical shaped energy density distribution for the laser beam 
as depicted in Figure 1. As shown, this configuration is most often used 
to size submicron particles. The mirrors are reversed for measurements of 
large particles. Thus, the open cavity region is bounded by a Brewster's 
window and either a plane or spherical mirror. Since the radius of curvature 
of the curved mirror is extremely large compared to the size of the beam, 
the geometry of the beam is nearly a perfect wedge with either plane or 
-spherical mirror (see Figure 2) . This wedge has a refractive index of 
1.0002922 for air or 1.000 for a perfect vacuum. It has the effect of dis- 
placing the beam on amount which can be calculated. The degree of misalign- 
ment that can be tollerated varies with the mode tuned. For TEMqq operation 
Bloom (1968)’'" gives a value of 5 minutes as to the angular allowance 
tollerable. For higher order modes, the margin increases approximately 
linearly with the degree of ordering. In our case, we must consider the fact 
that it is necessary to tune the laser at ambient air pressure prior to 
vacuum introduction. Thus, one is interested in the change in power output 
and the amount of mode shifting. 

The computed angular error due to a change of medium pressure from one 
atmosphere to a perfect vacuum is only 16 seconds. Thus, the only potential 
problem might conceivably be mode shifting. To examine the possible effects 

1 Bloom, Arnold: 1968, Gas Lasers . 
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FIGURE 1, Optical System Diagram for ASPS Systems 
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of mode shifting experiments were conducted at PMS to determine the amount of 

angular detuning necessary to produce mode shifting. A precision mirror 

gamble was used where differential micrometers could read displacements to 

10". The measured error required to shift from TEMq^ to TEMq-^ or TEMqq was 

p minutes and 3 minutes respectively. Thus, no problem would be 

expected from the effects of air pressure reduction to vacuum. 

PMS also analyzed data from vacuum ASPS operations at the Martin 

Marietta Corporation during an MSFC funded SKYLAB contamination ground test 

program (SCGTP) during May, 1972. Those results indicated that mode 

shifting did not occur at pressures reduced from approximately 0.8 atmospheres 
—9 

to 10 atmospheres. 
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IV. EFFECTS OF LASER HIGH VOLTAGES IN A VACUUM 

The high voltages required for laser ignition and continuous operation 
demand particular attention. It must first be remarked that the pressure 
range over which system operation is desired may encounter any of the basic 
ionization related electrical phenomenon generally classified as spark break- 
down, glow discharges (corona), or self-sustaining electrical arcs. 

First we will consider the nature of the expected laser operating 
potentials. Figure 3^ shows the voltage-current characteristics of the plasma 
tubes used in ASPS systems. Two characteristics are particularly worth noting. 
First, before the operating range can be reached a high voltage must be 
applied to initiate ionization. Second, in the operating region there is a 
negative resistance. The negative resistance region of the tube VOLTAGE- 
AMPERAGE characteristics require a limit to the supply current provided by 
~a*.ballas.t..,fresister -,of. .1-2-5K- ;ohms -*fco-.d50K -;ohms . The-'powar supply- block diagram 
of Figure 4_ shows a basic high voltage supply (1600V) with an add-on higher 
voltage doubler which has an added trippler network in series (High voltage 
ionizer network) prior to ignition and subsequent effective shorting through 
the ballast resister. 

Prior to laser ignition, the anode voltage nears 2500 V and this 
becomes the highest potential for consideration. However, it should be 
pointed out that electrical discharges often require a certain "charging" * 
time at the discharge potential prior to the actual discharge. This may be 
longer than the time required for laser ignition and thus preclude the external 
discharge. Thus, the actual ignition time is of some significance. To 
isolate the various aspects of the problem, we considered the three types of 
ionization phenomenon separately beginning with spark breakdown. 




LASER POWER SUPPLY BLOCK DIAGRAM 



LASER POWER SUPPLY CIRCUIT DIAGRAM FOR 1600 V. 


FIGURE h 
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As the electrical field between the anode of the laser and some 
external point (e.g,, nearest ground potential such as cathode or a mounting 
structure which is near ground potential) increases, the current increases 
steadily until some point there is a sudden transition from the TOWNSEND, or 
"dark" discharge to one of several forms of self sustaining discharge. This 
transition or spark consists of a sudden change in current across the gap. 

For plane parallel electrodes, the result is a spark that can initiate an 
arc discharge. For sharply curved electrodes, there may be corona or brush 
discharge. These latter discharges will be considered later. Under some con 
ditions the change in current may be quite small and the examples are typical 
but not singular in definition. For example, when a corona discharge 
occurs between a cylinder and coaxial wire this initial corona discharge, 
although a self sustaining discharge, represents an incomplete failure of the 
4 «gap ...ocGuring only ;in -a -limited., region -near .electrodes of small "radius of 
curvature, while the rest of the "gap" carries a "dark" current. In an 
engineering sense, breakdown is generally considered to occur only when the 
entire "gap" is bridged with the establishment of an arc discharge as common. 

1. Sparking-Potential. 

It was experimentally determined by PASHEN in 1889 that the 
sparking potential is a function of the product of the density and gap length 
only. However, most workers use the product of pressure and gap length. The 
sparking potential for air at one atm. pressure is given by 

V g = 30d + 1.35 (kV) . 

the spark-breakdown voltage has a minimum value at a critical value of pres- 
sure and distance (d) as shown in Figure 5* The important information here 
is the minimum sparking potential of about 325 Volts for a separation of 1mm 
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and pressure of 5.5 torr,or ^/2 torrand 1 cm. Since the likely separation 
is of the order of a centimeter or two, the greatest potential problem is 
likely to occur at a pressure between 10 ^ and 1 torr. 

2. Corona Discharge. 

The electric field strength (critical corona gradient) at 

which corona discharge begins is nearly the same as for spark discharge. 

In spite of the uncertainties involved in the derivation of the required 

field strengths for corona discharge, they compare most favorably with the 

emperical relationships given by PEEK (1925) . Corona is distinguished by 

the fact that increasing the separation between sharply curved surfaces, such 

as corners and points, generally results in the breakdown of the gas near 

the surfaces less than the spark-breakdown voltage. This breakdown 

appears as a glow discharge. While weak corona may not prohibit laser 

operat-i-on its persistence causes a .deterioration of insulating materials. 

Corona from an anode was generally characterized by bluish coloration over 

the entire anode while corona at the cathode (required grounding of the anode 

since the cathode normally runs within a hundred volts of ground) was 

characterized by red filaments of glowing gas along points and edges of the 

cathode. The critical corona gradient depends only on the radius of the 

anode (or surface features of the anode) and is not markedly affected by the 

effective "gap" distance. This independence of spacing is expected since 

for large separations the field at the surface of the anode would be relatively 

unaffected by further separation. This fact, that the critical field strength 

depends only on the anode radius facilitates the estimation of the minimum 

radius of curvature of metal parts used at high voltages. In the case of the 

ASPS lasers, the anode radius is 4 mm. The critical corona gradient is given 

for parallel wires by PEEK as 0.301 

E c = 30m b (1 + 8a ) (max. kV/cm) 



where b is the pressure in torr and T is the temperature in degrees Kelvin. 

It should also be remarked that the mechanisms of discharge about the 
anode is quite unlike that for the cathode. In the case of the cathode 
positive ions are formed and accelerate in the cathode field to produce the 
necessary electrons for a self-sustaining discharge. In the case of the 
anode, electrons formed in the surrounding space by various ionizing 
radiations always present accelerate in the direction of the anode establish- 
ing electron avalanches. These avalanches maintain the ionized state near 
the anode. 

3. Arc Discharge. 

The electric arc is a self sustaining voltage having a low 
voltage drop and capable of supporting high currents. This is the type of 
discharge already discussed with regard to laser operation. At atmospheric 
pressures the arc is characterized by a small brilliant core surrounded by a 
cooler region of flaming gases (Aureole) which is a region of chemical 
activity. Because of the high temperatures of the core, the gases are largely 
dissociated. At low pressures, the appearance of the arc depends upon the 
slope of the discharge tube. A constricted tube such as a laser capillary 
thus gives a bright glow even at low pressures. The absence of such a 
confining structure around the laser envelope makes it unlikely that arcing 
would occur external to the laser at low pressures. At higher pressures, 
the laser power supply simply cannot deliver the required current. Thus 
arcing was not considered in the establishment of the design criteria for 
laser operation in a vacuum. 

4. Summary of Design criteria and Final Design. 

The preceeding sections developed and discussed the theory 
relevent to potential problems with laser operations in the pressure range 
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from one atmosphere to a hard vacuum. The choices available to PMS for 

operating the laser were two: (1) operate the plasma tube in the vacuum 

and provide adequate sealed insulation around the laser electrodes to prohibit 

spark-discharge and corona, or (2) devise a seal for the laser plasma 

tube which would allow the laser electrodes to operate at ambient pressure. 

From tests it was determined that to operate the laser tube in the vacuum 

the electrodes would have to be potted in an epoxy or silastic material. 

(the minimum value in the Pashen curve of 325 Volts would provide for only 
1 cm separation at pressures of 0.1 to 1 torr) . Thus, one would estimate 
that a separation of centimeters would be required if the laser electrodes 
were not sealed over completely. A second point for consideration is the 
necessity to dissipate the heat as discussed in section 2. The amount of 
metal required to conduct away the heat places difficulty on establishing 
its safe separation from the electrodes. Thus, the sealed insulation of laser 
electrodes is again required so that this separation distance can be reduced 
without increasing electrical problems. 

One final factor that detracts from potted electrodes is the necessity 
to cable power from the power supply (at ambient pressure) to the laser 
(in vacuum). A connector vacuum penetration is required. If this connector 
is not totally sealed, its pins produce another potential spark-discharge or 
corona problem-*-. Such a connector if sealed would lose utility and re- 

a 

potting work would be necessary every time the laser had to be removed or 
replaced. 

In summary, the laser plasma tube was found to present a considerable 
problem if one was to insure its safe operation under vacuum conditions. On 

^Testing by PMS conclusively showed that Corona would develop 
between connector pins where laser voltages were applied. 
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the other hand, the change of medium from one atmosphere air pressure to a hard 
vacuum within the open cavity laser section would present no descernable change 
in laser operating characteristics. Thus, the logical conclusion to be 
followed was to devise a seal for the laser plasma tube which would allow the 
electrodes and bulk of the plasma tube to operate at ambient air pressure. 

Such a seal was designed and tested at PMS. 

For the seal to be safe, it was necessary to have special laser plasma 
tubes designed which had the cathode placed nearest the Brewster's window and 
with ample space between the cathode side-arm and the Brewster's window to 
epoxy the laser envelope into a stainless steel mount. This stainless steel 
mount was constructed with a VITON "0" ring seal which is compressed using 4 
socket head screws. The epoxy can be most any commercial clear apoxy having 
a long curing time. This long curing time allows for complete out-gasing and 
.the best laser to metal seal. .It is the opinion of PMS .that this seal is 
superior to the standard seals used in the construction of the laser itself 
and thus is not expected as a point of failure. 
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PART III 

Theoretical Response of Active Scattering Particle Spectrometers 

The scattering of electromagnetic radiation within an open cavity 
laser was recently modeled by R. G. Knollenberg at PMS . A computer program 
was devised to compute the scattered energy for spheres having variable radii 
and refractive indices in a plane wave. It was recognized that in the ASPS 
the collected light signal consists of both forward and backward scattered 
light and that these can constructively or destructively interfere. However, 
after further analysis it was recognized that the only time the backward 
scattered light accounts for even 10% of the total light signal is when the 
particle size is considerably less than the wavelength and under such con- 
ditions the phase difference between the forward and backward scattered 
light is nearly the same. Thus, the MIE scattering computations are 
adequate in their classical form and the necessity to solve the problem of 
electromagnetic scattering in a standing wave was precluded.^ - 

The computation involves only four input par am eters -wavelength (X) 
radius (r) , refractive index (m) and scattering angle (0). The particle 
radius and wavelength are generally combined into an parameter x given by: 
x = 2 ir r/X 

* 

Because of the large variation in scattered energy (intensity) with angle, a 
logorithmic plot is used. A sample of such a calculation is shown in Figure 
1 for a particle radius of 1.0 micron and a refractive index of 1.6 (Real) 
with no absorbtion (imaginary part of refractive index = 0. 

^To this writer's knowledge, the problem of electromagnitic scattering 
in a standing wave has not been solved theoretically in exact form. The 
boundary conditions would appear not identical to those given by MIE. 
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In addition to this , computation, the total scattered power was 
integrated for each degree step to enable the determination of the amount of 
scattered energy received by the collecting optics in ASPS systems. Printouts 
of these incremental scattered intensities are shown in Figure 2 for scatter- 
ing angles of 1-28 degrees and 152 to 180 degrees. The quantity ITSUM is the 
total scattered energy collected over 4 tj steradians. It serves as a normali- 
zation parameter and can be used to estimate the efficiency of the collecting 
optics . 

Note the uniformity of scattering signal vs. angle over collecting 
angles up to 20° ; The sum over the collecting aperture of 5-20° assures a 
smooth calibration for ASPS instruments at least in the vicinity of one 
micron as shown here (see figures 3a & 3b). At radii of 0.1 and less, the 
refractive index can cause a change in scattering power by nearly 4. Choosing 
a refractive index value of 1...6 allox-js for an uncertainty of ±2 in scattered 
intensity for completely unknown refractive index. This results in an error 
of ± v 2 in terms of particle radius which would account for much less than 
the specified rms error of ±0.1 microns. If the refractive index of sized 
particles is known, one should use the theoretical computation for that 
refractive index for the most exact results. 

We have run computations for real refractive indices from 1.3 to 2.5 
at intervals of 0.1. The computations are made every 0.1° and for radius 
intervals of 0.1 micron for radii up to 100 microns. The resulting computa- 
tions are filed on 90,000 frames of microfilm at PMS and are available from 
PMS. The entire program is listed on the following pages. This program 
requires a considerable amount of machine time and certain plotting routines 
are not generally available at most institutional computation centers. 
Questions regarding its use should be directed to R. G. Knollenberg. 
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Computer Printout of Energy Scattered into Incremental 
Solid Angles, Notes ITSUM = The total scattered energy 
collected over 411 steradians, 
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Computations as just described were used to calibrate the ASPS size 
ranges. This theoretical calibration is approximate in that it does not 
account for particle shapes other than spheres. The selection of the collec 
ting angles minimizes the effects of irregular shaped particles. In the 
large particle size range, calibrations were performed with polystyrene and 
glass spheres as well as samples of irregular shaped particles. These all 
support the theoretically derived calibration curves. 
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4 I.M A-G-I-N AR-Y PART E __ * F. 1 Q. • 6 _ *___ J_* _ /.J 

85 FORMAT! 3X , 16HSCATTER IMG ANGLE * 1 OX *2H I 1 * 15X t 2H 1 2 » 15X » 2HI J ♦ 15X t 

4_XWLAjXAX_>AIlL6yXb9X_OXHPOLARIZiX_ION__/__9X_XlIt , §JAl__ / J T 

95 FORMAT ( 8X »F6 .2 » 5X ♦ 6E17.6) 

105 -F LDRMAJ 1.A6.X-, *S 1 7 F PARAMETER = * F1Q.A / ) 

115 FORMAT ( 1 0 X *N* 7X 2(7X *A SUB N * ) 3X 2(7X *3 SUB N*) 9X * ( A SUB N 

44X.Xb*JL7_6.X_XLB_XUX^NLLl_3i± > l*2_XJXXXL10Xj^?P.^.t_§^.__lIPA G J^. A 3_ Y _‘ ! ll. / !_ ) _ 

125 FORMAT (9X » 13 »9X «1H( * E11.4,3H » ♦ E 1 1 • A ♦ 1H ) ♦ AX ♦ 1H ( »’E 1 1 .4 * 

1 3H > * EJU_.A.*mLiIX» Fi_L._Ajl_X > fJLLifjJ 

135 FORMAT (*-* / ^-EFFICIENCY FACTOR FOR EXTINCTION - * E12»&) 

format t » — FFFIC1ENCY FACTOR FOR SCATTERING = * E12 * 6> 

155 FORMAT (* — * / * — * / * — * 10X *TFI1S PROGRAM WAS PROVIDED THROUGH The 

j CO URTESY OF .PARTICLE MEASURING SYSTEMS* BOULDER » COLORADO* ) 

165 FORMAT (#1* 35X *A SUB N AND B SUB N COEFFICIENTS OF MIE FUNCTIONS 

.1*-/-) : 

175 FORMAT ( 3 X ^SCATTERING ANGLE* 7X ^AMPLITUDE S SUB 1* 1AX *AMPL I T UD 

lj f__S_S UR — QX.*T H FTA* ?X ?(..] 0X * REAL * 8X.,*lM AGI NARY* ) J _ ) 

185 FORMAT ( 8X »F6 .2 »7X *1H( » EU.A*3H , » E 1 1 . A » 1H ) » AX » 1H ( *E 1 1 . A » 

1 — 3H — » — *. EL1-L*AjJJ±L) 

195 FORMAT (*1* AOX ^COMPUTATION OF MIE SCATTERING FUNCTIONS* /) 

__ 2-0 5— F-0 RiVLAJ. L6 E.lil. 2., t.i 4_,_L4 * AJ± _M J EJ 

215 FORMAT ( AE15.6 *8X , I A ? I A .AH MIE ) 

_22-5--E.QRMA.T-X- ... 3-EUL»Aju I A X v I .A. A XjlA±LM JLSJ 

2 a 5 FORMAT (*1* 2 OX *AMPL I TUDE FUNCTIONS FOR MIE SCATTERING* /) 

2 6.5 F.0 Rf iA-L— l 1X-*J_FEDE X_PF RE_F_R AC T_ION 1_J5 JlAL_ PAR L = 

1 IMAGINARY PART = * FI 0.6* * >* /> 

275 FORMAT- ( 15X <-SIZ£ PARAMETER = » F1C«A, 1UX »DTHE T = » F5«2 /) 

285 FORMAT {*-* / -“--* / *- ELAPSED TIME FOR THIS CASE = * E10.2* 

1 * SFF ON OS . » ) 

310 FORMAT ^-NORMALIZATION FACTOR = 

3JL5.„EO-RiM_A.T (_*!*.) 


2 


015 


02 r " 


626 ” 


035' 


038 

”039 

0A0 


* L12.6) 


WRITE(6*155 ) 

RF AD ( 5 . 5-Q.0J--L.Ll3 ULX a L.T Y_»L I T * L I T T 

) FORMAT ( 4 A 1 0 ) 

“358 97 932384626A 3D0 


CONST 2 = CONST1/180.0DU 

D _TJT ET..I =JL« - 

XLAM=.6328E-6 

__READa5..,Xai_Pil,DjHEJ_A J NJJlF_T 1 RX 

) TIME IN = TIMEF(O) 

U-SUFIXj 

RX=RX+0. IE-6 

IJEL _( R X.. £ Q_*ill_Q O j E_-6 i_ G_0 _ I P_ J,0 9 2 

X = 2 -'CONST 1*RX/XLAM 

E.N.C OD.EJ. 1 .0_»_5_Q. 1 j_L_LI 1 4 J J X 

ENCODE < 10*501 *L ITT (A ) ) X 


WRITE! 6*203 ) X 

3_£-0.RMAJXZ.2./.»3_Q.X..iT-X = !_*_F_LQ-»J3J 

DTHET = DTHETA 

.I.F (-NJH E.T.1- Jj03.a1 0_0_0 

3 IF (NTHET .GT. 181) GO TO 1001 

R = REALJP_iil_ 

RP=AIMAG(PM) 

LEM P.RP _ =.. JR P 

R P = -R P 

Q_T d =_I .*..2. *.X±.9_ «. .0. 

SN = S I N ( X ) 



A IMA = HSIN*HC0S/TEMP3 

AO-=C MR LX. { ARE A LXU M AJ 

TEMPI = 1.0/X 

_ _WX -lJ— = _ J-EME l^WOrWM 1 

W ( 2 ) = 3 .0*TEMP1*W ( 1 ) -W'O 

-2LG-MR2-C — s— LEL1R-1XEM _ 

A ( 1 ) = - TEMP 2 C + 1.0/(TEMP2C - AO) 

__W R EAL.= BEAL.I.WJJLLJ. 

V/RE ALO = RE AL ( WO ) 

.. TEMP3C . =.. A ( 1 ) / PM +_ JXMPJL 

TEMP AC = A ( 1 ) *PM + TEMPI 

SAi-1-1 s i TFMP^C*WREAL - WREALO ) / < TEMP3C*W ( 1 ) 
SB C 1 ) = ( TEMPAC*WREAL - WREALO )/( TEKPaC*W l 1 ) 


DO 102 N= 1 *400 
N2=N+2 


PN 1 =N 1 

XEM 9-2— s— (-2 *N 2- 1 J JtT£.MP-l 

W(N2) = TEMP 2 *W ( N 1 ) - W(N) 

— W4RE-A-L-^R&AL-. 1 -W-UZLU 

WREAL 1 =REAL ( W ( N ) ) 

_T.fE-MP.3C =--B N-1-&-T-EMP-2C. 

A ( N 1 ) = -TEMP3C + 1 *0/ (TEMP3C - AIN)) 

-T.E-M P 3— =~ EM -1JLT-EM.PJ 

TEMP AC = A ( N 1 ) / P M + TEMP 3 

T E m P^C-5— AXN-1J-£E M— ± J.E.MP.3 

SA (-N ] ) = ( T E M P 4 C W REAL - WREAL 1 ) / ( TEMP4C*W(N1 ) "WIN ) ) 

SB-UM-1-) ^---LT-EMiP.5_C.*w.RE^L_r'_ZtBE Ai^Ll/JXEMPj'iL^iilLL-.ri^J^.ii- 

RSA ( N ) =REAL ( SA ( N ) ) 

; R-S B-( -N-)-S RE. A L-L S.B-LN. X.) 

URSA (N ) =A I MAG ( SA (N ) ) 

UP-Sfm\U--AJJiiAjGuLSBJ-NJ-J 

ACON J = CON JG (SA ( N ) ) 

3-GO.N J=jGQIsI J.GJ-S.9. (El J_L 

SA2(N)=SA(N)*AC0NJ 

S R-2- C-N-)-=-SB .t C Q.N. J. 

RSA2 (N) =REAL (SA2 (N) ) 

R^p.?-UPP^47-AB4-SB.2-LR-U 

TEST 1 = (RSA2 ( N ) + R S B 2 (N) ) /PN 

I-P4-TE-S-T-1---T-ES.LZ)-3. J 2_»-Z 

2 IF (PN-OTH) 102*3*3 

- 1- 0-2- -CO N T-I-N-U £ 

WRITE (6 *55) 

3-KBBRN-s.N 


-9_9__D0 —5 — I — - — 1-*N IHE-T. 

THE T A ( I ) = FLO ATE ( 1-1 ) *D THETA + DTHETI 

T-H i-T-(-I-)— =— T-H E-TJU- 1.) 

IF (THET ( I ) .LE.180. ) GO TO 7 

NXUSXaliJ 

GO TO 6 

J CT.dEJ-A.--=_-STlGL.(-D_.C-0._SJ J_HEI Aj_Lli COJiST 2_)_) 

C THTS ( I ) =CTHET A 

S.T. HE J-A. - =. . ISN.G L i .D.S.I N J J.H£ I AJLL K CON S 1 2 JJ. 

SKI) = CMP LX ( 0. *0. ) 

9 2 i I- ) =lC MPLX_LQ-.-.-Q-«Q 

PI (1 ) = 1 - 

R .U.2 )_ = 3-« -O-f.C-T ±t£-L A, 

TAU ( 1 ) =CTHFT A 

T. EMP_4.---.-S-IHE.TA** 2 

T AU ( 2 ) = 3.0*(CTHETA**2 - TEMP4 ) 


DO 4 N= 1 * KEEPN 

_L'U_=N.+.l 

N 2 = N + 2 

...T E MP-2 — =— 2 * NZ--LL 

TEMP 3 = N 2-1 

pun?) = ( TFMP?*P I < NM *CTHETA - N 2 *P I ( N ) ) / T EMP 3 

T AU ( N2 ) = CTHETA* IP I ( N2 )-P 1 ( N ) ) - TL.MP2*T EMP4*P 1 IN 1 ) 
_ -T-EMP 2-=-2MtI 


ID. 


TEMP 3 = N*N1 

TEMP ° - tf.mp?/tfmP3 

5 1 ( I ) = TEMP2*(SA (N )*P I (N)+SB(N )*TAU(N ) )+Sl ( I ) 

<;2(T1 = TEMP 2*J_iiB. 1 NJ *P.I (N) +SA ( N ) *T AU ( N ) )+S2( I ) 


4 

3 

CONTINUE 
CONI I-N.U-E 

142 

C 

6 

pip 12 1 = 1 .NTHFT 




S1 RE( I )=REAL(S1( I ) > 

51 IM ( n - A I MAG ( .91 (111 __ 

19C 

191 



S2RE< I ) =REAL ( S2 ( I ) ) 
32 IM ( I ) -A I-MAG ( .3?„.( I ) ) 

192 


12 

CONTINUF ^ 


£~ 


COEPH = 2.0*TEMP1**2 

DO 8 N = 1 » K.EEPN 
TFMP? - ( 2 * N + 1 ) *C0EPH 

— 

EXKROS = TEMP2* (RSA(N) + RSB ( N ) ) + EXKROS 
<;ricpn^ - TEMP o » ( RSA 7 ( N ) + R5B2 ( N ) ) + SCKROS 


8 

CONTINUE 

AN OPM ~ 4 . / (SCKRQS*X**? ) . _ 


C 


■ DO -9 I- - 1-.NIHEX . .. 

S1C0NJ=C0NJG (SI ( I ) ) 

FNT 1 ( U -S-U 1 1 * 3 1 CON J 

' 164 
162 



S2CON J = C.ON JG < S2 II)) 

EN T 2 1 I ) = 92 ( T )*S?CONJ . „ .. . . 

16c 

167 



RENT ! ( I-) =-REAL ( ENT HI)) 
PFNT9I i )=RfaL(FNT2< I ) ) 

168 

16S 



ENT3U )=REAL(S1( I )*S2C0NJ) 

FMA 1 T )=-A IMAGES! ( I )*52C0NJ) ' 

1 7 C 
171 

TEMP 2 = RENT 1 (I) + REN T 2 ( I ) 

pm a p i T ) = (RENT 1 ( I 5 - REN T2 ( I ) ) / T EMP2 

ENTAVMI) = TEMP 2/ 2.0 
r,0 TO 9 

IT ( I -1 1 = ( ENTAVG ( I ) +ENT AVG ( I-I ) ) *CONS T 1 * ( C TH TS ( I - 1 ) — CTH T S ( I ) ) 
jF ( I T.( T--M .1 F. 1 *F-1 2 1 I T ( I -1 ) =1... E- 1 0 

I B I S - I - 1 

_ WP I T.F 16_._2_D.2_) I B I S.».I.T L I.-_L ) 


202 

FORMAT ( 1H0 » 1 4 *E15 .6 ) 

I T.5 U M=I J.S.U.M -+-1 T LI.-1J _. _ . 



9 

CONTINUE 

ENCODE ( 2-0 5.9000 5 1 .1 T ( 3. 1 ) IT SUM 



5 000 

FORMAT ( * I SUM = *E 1 5 • 2 ) 

JJMF-OJ - I.LMEE1.0J 

DELTI.M = (TIMEOT - TlMEIN) / 1000. 0 
Ga T.O 9.9.99 

— 

c, 

c 


NO. t ni IT PI IT 


c 

— 

WLR.LT E-( - 6 - 1 - 1 6.5 ).. 

WRITE16.73) R*TEMPRP 

-WR ITE ( 6 .103 )X. . ... 

2 C. 



WRITE ( 6 » 1 15 ) 

WP 1 T p ( 6 . 1 v F ) (N.RSA(N) » URSA ( N ) »RSB(N) »URSB(N) » RSA2 l N ) *RSB2 ( N ) . N = 

21' 


1 l.KEEPN) 


— C- 




NO. II OUTPUT 


WR ITE<6 *245 ) 

-W-R4 XEX6 -,-2-6-3. )_ R-,XE MP.RE ; 

WRITE16»275) X. DTHET 

-WR4XF-(-6-*47_54 

WR ITE(6 *185 ) (THET ( I ) .SIRE ( I ) *S1 IM( I ) »S2RE( I ) *S2 IM( I } . 


1 .NTHEf ) 


211 

22C 


-G 

C NO. Ill OUTPUT 

-€- 

WRITF16.195) 

_W-R4XEX6-.,7-XL-R-,XEMP.RP- 

W R I T E ( 6 , 6 5 ) X,DTHET 


WR4-T -g-F4-»-3£4 : : — — - 

WR I T E ( 6 » 95 ) (THET ( I ) , RENT1 ( I ) , RENT 2 ( I ) »ENT3( I ) >ENT4 ( I ) » ENTAVG ( 1 

j.._’ 4?-Q.L-A'R-(-I-E. — I — =— -l-s-NXHEJ-) 


) * 


9999 CONT I NUE 

<:.AL-U- IOlO.T-4J-Hei3 ENXAyjGL*.N-T.hiEJ-*.2jiJJJLI8jJLj_X ♦ LJ_V iLJilliJ 

CALL DASHLN( 14308 ) 

CALX-CURV-E-l T.H E-T. t-REN T 1 i M T HEJ.) 

CALL DASHLN ( 1634B ) 

GALL. CURVE_(-Id.EJ_»-R EEEL2 aN.TJ8.E_I.). 


NTHET 1 =NTHFT -1 

CALL-4 0- 1 OJ4-IH EX,XX.vNXd EXL»_2 j . 1 JJ_7 8_, L.T X .LJ X t L.I J_T_._- 1J 

DTHET I=THETA (NTHET ! 

}-F4-T-H&J4-NJ-lrlEX U-I-I-.4.&0 .-)■■ -G.Q— IQ_ £2 

WRITE (6, 50 11) ITSUM 

§ 044- _ fo rm A J-( -* 0. - _ -I.-T-SU M=_ _* EX5_.,6 .). 

. NXHE-T = 1 8 IS DTHET 1=0. 

. GO — T-Q — 1.Q-0. 

WR I T E ( 6 » 1 3 5 ) EXKROS 

W-R44-E4-6-.4-A-54SCK-R.QO 

WRITE ( 6 » 310) ANORM 

WR4-T-E — ( -6-*— 2-83- ).— D.EXX 1 hi 

GO TO 100 


C- 

C END OF RUN 

G 

1000 WRITE (6 .3 15) 

GO— 1 TO— 100-2- 

1001 WRITE (6. 40) NTHET 

— , -1002— CON-T-I-NUE 

401 0 CALL FRAME 

— CALL — E-XJ-I 

END „ , r . 

SliRROilT INF IDIOT (FLDX .FLDY ,N ,LTYP£ , IPATRN ,LABX »LABY ,LT I T ,LFRAjvX ) 

C SUB ROUT I NE IDIOT ( FLDX » F LDY ,N»LTYPE»IP ATKN » LABX » LABY , LT IT »LF RAME , M ) 

C -SUB.ROUXl.NE 1 D IPT ( FLDX .FLDY , L , L T Y P E > IP AT KN,LABX,L A BY 4 T I T , LF RA Me , N ) 

C FLDX FIELD OF INDEPENDENT VARIABLE 

G FX-&Y £JXXD_QF— PEP F. NXE NX_Al ARJLAi3 L E 

C N NUMBER OF DATA POINTS 

G -M N UM 8.E.R _ -O.F C _U R.V.ES 

C LTYPE 1 LINEAR-LINEAR 2 LINEAR-LOG 3 LOG-LINEAR 4 LOG-LOG 

C LP-A T RM E J EL D _ OF DA S H PATJ.ERNS.. FOR. CURVES 

C LABY FIELD OF ALPHA CHARACTERS FOR FLDY LABEL 

r |_t;t FI FID OF ALPHA CHARACTERS FOR TITLE 

C NDX IF 0 FLDX IS ONE DIMENSIONAL IF 1 FLDX IS 2 DIMENSIONAL 

C I.HGP.S -X -GRID I.S_ HALF AX = 0__GRID 1_S _ G RUX _ _ J^l _ G R I D I _S_ PE R I M L 


187 


197 

19fc 


2 32 


. 12 . 


C I LAB = 0 NO LINE LABELS =1 LINE LABELS A = FLDY ( 1 ) B = 2 ETC 

,j^44GN-S4QN— &LLXX44La4-jll-0.-YlJLs.aJ.j.i^.J.l^^ 

C DIMENSION FLDX ( L » 1 ) .FLDY(L*1> ♦ I P A TRN ( 1 ) ,LABX( 1 ) * L ABY ( 1 ) »LT I T ( 1 ) 

D-I^E-NS-ION._I.R-E.C(.42j-A.Q1.vLS.Y.M.(Z6J...jJ_Q.L26J 

DATA ( LSYM=] HA . 1HB * 1HC » 1HD * 1HE * IMF * 1HG ♦ 1HH * 1HI * 1HJ * 1HK 1 1HL » 1HM » 1HN * 
J.H_W il H Y_»_lHZj 

I HGP =0 


-C 

NDX = 0 
M" 7 






c 

M= 1 

-N SF-A R-C H-— N- 
NSEARCH=N+2 

in k' -J . 





— 


20 T RFC ( K ) =0 

vw i w-YMAy = pi ny M .1 > ... . .... _ . 


YM IN=YMAX=FLDY (1*1) 

-DO- l -4=-l- »M_ 

DO 1 I=1,NSEARCH 

VM T W -M T Ml F ( YM 1 N . F 1 D Y ( I , 

J ) ) 

— 

— 


— 


ymax =max 

4F ( NDX. E 

IF ( YMAX i FLOY l I » J ) > 
Q.e.0-.AND.. J . G.T. . .1 ) G.Q J.0 

1 

XM I N = M I N 
X M A X- M AX 

IF ( XMIN .FLDX ( I * J ) ) 
1-F ( XMAX ., FLDX. ( I » J.Ll 


1 CONTINUE 

t f t op RAM 

F . L T , 0 ) C A 1 1 ...FRAM f 


IF (LTYPE 

. EO • 0 ) GO TO 8 
4*5.6) LTYPF . . 


-3 CALL L TN 
CALL 1 I N 

RD ( X M I N * X M A X * X B ♦ X T * MG 
RD ( Yd 1 N_s _Y jM AX_*_Y B >_Y T ,_MG. 

RX » 1 NRX vl FMT X »NUM-X * I XOR ) 
RY , INRY , IFMTY ,NUMY * I YOR ) 

GO TO 7 



4— GALL ■ L4-N 

RD (.XM I N-. XMAX *XB *JXT_>_MG. 

RX , I NRX * I FMT X * NUMX > I XOR ) 


CALL LGRD ( YM IN * YMAX » YB » Y T »MGRY » I NRY * I F M T Y »NUMY » I Y0R ) 

=.60—10 — 7- 

5 CALL LGRD ( XM I N ♦ XMAX *XB *X T *MGRX ♦ I NRX » I F MT X *NUMX ♦ I XOR ) 

.CAL- u - .LAN R D-LY.M I tL*_Yj- i A X j_Y_B. .Y J_*.M G RY_.*_ I N RY_>_ }_F MJX MY ] Y OR. [ 

GO TO 7 

6^-GAA A-4XAi0AJ0MAN-»-X^ L.I L MI X »NUMX UXOR. }.., 

call lgrdi ymin . ymax ,yb .yt ,mgry . inry. ifmty .numy , I YOR ) 

7-_-C-/vLL--SFJ-T-(«-1-5-»-. .9A-»_..15.j-«-9.0j.XB_>.XIj„Y.B.a YJ_»_L IYP_E_L 

CALL LABMOD ( I F.-iTX ♦ IFMTY »NUMX * NO MY *1»1*G»U»IX0R) 

S — I-F-LI-HG P_).9..*-1AL».L1 

9 CALL HALF AX ( MGRX ♦ I NRX *MGRY * I NRY »XB *YB > 1 » 1 ) 

GO— T-O— 1-2 

10 CALL GR I DL ( MGRX » I NRX »MGRY » INRY ) 

GO— T-0 12 

11 CALL PERIML (MGRX » INKX.MGRY. INRY) 

1-2— MM =4-00-4 XOR *S.a 

CALL PWRT ( 24 5 .MM.LABX >40*2*0) 



CALL PWRT (MM, 220 ,LARY *40,2 , 1 ) 

GAL L— P-W R-T-L60 ., 5 6-5-*. L J_I.T.-,.4.a_ J _3.,D J 


JJ = 1 

DO — 1-7— J— 5- »M 

IF (NDX ) JJ = J 



13 DO 14 I = 1 » N 

14. _CAL! PLAN I .( F.LD.X .( .T.r.J.J. ) _».F_LD_YJ J_»_J .) J 


GO TO 17 

-*5— eyVtrtrH^rS^lrN-<^^^-M-0-M 

IF ( IPATRN ( J ) .LE . 1 ) CALL DASHLN ( 1 7 7 7B ) 

(C-AT^L— FR^S-T-FT-I-F-L-OX-I-L^J J4--> fLLD-Y-(J_.»-JJJ 

DO 18 I--1>N 

GAL-L— VEC-T-OR-(-F-L-DX-(-LvJJ-)-*-F-LD.YJ_L>.JJJ 

CALL MXMY ( MX »MY ) 

1X= < - MX — 130 - F/ - 2 0 

I Y = ( MY- 130 1/20 

^,f_(_I-X^.L-T-*-l-*OR-*-I-X-»-G-T-»Al-.-OR»-IY-»-L.L»l_«_QE.».i_Y_»J3I — GO — IQ — 1_8 — 

I F ( IREC ( IX , I Y ) ) 33 ,31 ,32 

-3-1—1 R-E-G (~I-X-,4-Y~> -=-J 

GO TO 18 

-35 — Ff-H-R-FG-F- TX G.Q - . I O., . 1 8 

3 3 I F ( I R E C ( X X * -I Y ) ,».L T • - 1 0 0-0 00000) G O T-0.18 

FR.g G-I-I-X- ■* 4 -Y-)-=- =-( 4-0-QS X A ELSE A _LR.EC 1 _UC ». 1 Y..LL+ J J 

18 CONTINUE 

- F-7 — C-G N T-I N U E- 

CALL DASHLN ( 1777R ) 

-FF-UM-o-E-G -o-l-) — 00 — T-0 — 8-1 

IF ( I LAB ) 80*81 

--8-Q-- GON-T-I-N-UE- 

DO 5 0 K = 1 , M 

-I G-LK-l-= 0 

DO 60 1=1,40 

-90— T-O-G-^l-^3-8 — 

IF ( IREC ( I » J I „£Q.K )65 ,70 

-7- 0- -GO N-T -I-N U E r 

GO TO 60 

,_6§— GALL— R-SYM-(-lA.0^2.0*-i-»-lA.Q.±ZD_* J.aL_SXM11>J_».1 

IG(K)=IG(K)+1 

HM-+-1-0 

6G CONTINUE 

— § 0— GON-T-I-N UE 

81 CONTINUE 

]-F- (• LF-R-A-M g-.-G-T- *0-)-CAL I — FJ3. A M E 

RETURN 

SUBROUT INE LINRD(XM1N,XMAX,XB»XT ,MGRX , INRX, IFMT »NUM» IOR) 

4-GN L 0 G— - ..4-342-94 A8-L-9.03.3 

RANGF=XMAX-XMIN 

I . F L(-RANG&.L^AP5F-LXMJTL*L,F-_ai-JEAj^G£=.XidLN*.iIl 

TF (RANGE ) 10 , 1 1 * 1 0 


X T = 1 . 

MGRX-— 1 — 

INRX =2 

I-F-.VJ—6-HT -E-4-..1-) 

NUM = 4 

FOR-=-0 

RETURN 

- 4-0 — CO N-T-I-MU E 

I LGRG=TENLOG*LOGF (RANGE ) 

I F (RANGE . LE . -l-.-U LC--RG=J.LGEGrl__ 

R0NDRG=10.**ILGRG 

R9 N =94 (NGE-ARCLiLS-G 

RDX = 1 • 

-1 F-l-R D N . LE- . 6-.-I-RD X=-«5. 


IF ( RDN.LE . 2 . )RDX = .2 

B+V-54-Z-^R ON-DRO^S-D* 

J1=XMIN/DIVSIZ 

I-F-(-X M-I-N-. L -T-.-0- )- J4-=-J-L=-l 

J2=X MAX/D I V S I Z 

jf.p_(_X ,v) A X-.-G -T-.-0 - )- J 2-^ J-2.+ 1 , 

IF (J2.EQ.J1 > J2=J2+1 

XJ3==0-1-JUAMV-G4-Z ; 

XT = J2 *D 1 V S I Z 

MG-R-X—-J-2---J-1 

I N RX =M I N 1 F ( 10.*RDX+.1 *5.05 ) 

— —}.p_ (- J. L-GR^v .. LF-.-Q4 -I. -_LL.G RG_+_1 

IF ( ILGRG.GT.U) I 0 = 0 

7^-04-V^-I^Z — 

I F ( ABSF ( XB ) .LE . EPS ) 1 *2 

— 1 — X-R-=0-» 

M I N L = 2 

— - — G Q — T 0 — 3 

2 M I NL = TENL0G*L0GF (ABSF ( X B > ) +1 

LI . ..l - 4-M j-&L=4 

3 IF (ABSF (XT ) ,LE.EPS)4»S 

— 4— X-T--0-. 

M A X L = 2 

GO— TO — 6 

5 MAXL=TENLOG*LOGF(ABSF(XT ) ) +1 

pf^(J4-A-X-U-»J=X^l-)— M aX-Ls-1 

6 IF (XB.LT.O }MINL=MINL + 1 

1- FJ-X-T- .L-T-.-O.AM A X.L=MAX.L ±1 

NUM= I D+MAXOF (MI ML >MAXL ) 

1 F4_M5.;.N JJM-=.N UM+-1- 

IF (NUM.GT, 10)8*7 

— 7__p;4GG)(^!-7-.4Z>G^4X4Vaom i :Z-i-LD 

00 FORMAT ( 2H( F » 12 * l'H. » 1 1 »1H ) ) 

GO— T-0--9 

8 IF ( ID.GT.3 )ID=3 

MUM — — 7-+-I-D 

ENCODE ( 7 , 101 » IFMT )NUM * ID 

ra/| A44-2444^.4^^1^ j J_U-1±LU 

9 ICRT=12#NUM*MGRX 

1-0 R-= 0- ^ 

IF ( ICRT .GT .790 ) I0R=1 

RF-TUR-N 

END 

GUBR-OO-T-l-iLF L-ORDJ JGZ.I-ui-.X-iL AXj XB_r XJ. .jJilG-RX jJLNLK X. » i_FMLjj X U. illJ.Qiil 

XM I N =XM I N * 1 . E- 1 S X M A X = X M A X * 1 • E + 1 

TE-N LC G-=-«-4.3 4.2-94-4.8-1-9-0.3.2 

I XB = TEMLOG "LOGF (XMIM ) +.001 

I-X-T— T-E N LO-G * L.QG£-( XMAX _1__t . QDJ 

IF ( X M I N . L T , 1 • ) I X B = I X B - 1 

PF-I-X^MA-X-. GE -^ I .^- U - X -I =_LX.LJJ 

IF ( IXB.EQ. I XT ) IXT=IXB+1 

XR--lrO*^-*-I-XB 

XT=10.**IXT 

MG R X = ( I X T - 4 X 34- / 2-C +-L 

I N R X = ( I XT-1 XB ! /12* J o. 

MUM-=£ 

I F M T = 6 H ( E 8 . 1 ) 

I-G RT-=-(-J-X J -4-X.B J *9.6 


IF ( ICRT.GT.790) J0R=1 

-ftfrWfcN 

END 


— -K-R-UN 

IAVG (LOG) VS THETA $S X = 

-T-H6T-A— <-9 EGTTE-&-SH 

AVERAGE INTENSITY (LOG) 

EGON — (-LOG-) 

ICON VERSUS THETA 

-0-.-eQ--l-.D-C-- - 1 - 8 - 1 - Q.Q...E--^6- 

0 

1 -.- 4 Q _a.-GQ--l-.-DO 1 - 8.1 aa..£_- 6 ._ 

1.60 0.00 1 * DO 181 00.E-6 

i-^7-0 .. 0 . GO - l— -D-C 1-8-1 Oa. iL= i»- 

1.80 0.00 1 . DO 181 CO.E -6 

. 3 ... G-.-0-a-i-.-D0 — i- 8 -i 0-0. . e^-js,. 

2.00 0.00 l.DO 181 00.E-6 

2 -.- 1-0 0-.-Oa--l-.-DO -1-8-1- -0-Q-.E--6-. 

2.20 0.00 l.DO 181 00. £-6 

- 2 -.- 3 G, O-rOD— 1 -.D 0 1 - 8-1 aa.£-=i>- 

2.40 0.00 l.DO 181 OO.E -6 

0-.-00-- 1— DO 1-8-1 0-0-..E-— 6 _ 


*END 



